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ABSTRACT 


Hypervelocity  flows  for  velocities  in  excess  of  1.4  km/sec  (Mach  5)  require 
very  high  stagnation  temperature  to  avoid  liquefaction.  The  arc  heater  wind  tunnel 
has  been  designed  to  provide  such  flows.  The  electric-arc  driven  wind  tunnel  can 
develop  stagnation  temperatures  up  to  13,000*K  which  will  produce  hypervelocity 
flows  up  to  7  km/sec  (earth  orbital  speed).  The  nature  of  the  flow,  however,  is  such 
that  the  high  temperature  source  flow  may  cause  severe  gradients  at  the  nozzle  exit. 
In  order  to  perform  aero-thermodynamic  tests  the  characterization  of  the  flow  in  the 
test  section  is  required.  This  paper  experimentally  determines  the  stream  profiles  for 
an  arc  jet  wind  tunnel  conical  nozzle  directly  from  calorimetry  and  pitot  probe 
surveys. 
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I.  INTRODUCTION 


This  report  is  concerned  with  developing  a  procedure  by  which  flow  field 
characteristics  of  very  high  temperature  arcjet  flows  may  be  produced  directly  from 
stream  surveys  of  nozzles. 

Hypervelocity  flight  in  the  atmosphere  has  previously  been  concerned  with 
entry  into  the  atmosphere  and  descent  to  the  ground.  This  has  been  expanded 
during  the  last  decade  to  include  vehicles  which  will  use  the  upper  layers  of  the 
atmosphere  to  perform  hypervelocity  aerodynamic  maneuvers  and  vehicles  which 
will  fly  from  the  ground  into  orbit  using  air  breathing  powered  lift.  The  design  and 
testing  of  an  aeroassisted  orbital  transfer  vehicle  (AOTV)  and  the  National  Aerospace 
Plane  (NASP)  will  require  ground  based  testing  in  order  to  obtain  a  data  base  for 
this  flight  regime.  Such  data  will  be  used  to  validate  3-D,  reacting  real-gas,  CFD 
codes  that  will  be  used  in  the  design  of  advanced  vehicles  where  aero- 
thermodynamic  effects  are  important. 

The  60-MW  Shuttle  Interaction  Heating  Facility  (IHF)  was  designed  to  meet  the 
testing  requirements  for  development  of  the  thermal  protection  system  of  the  shuttle 
orbiter  vehicle.  Its  design  provided  for  heat  transfer  simulation  and  stagnation- 
region  flow  simulation  for  conditions  which  the  shuttle  orbiter  would  encounter.  If 
this  facility  is  to  be  used  in  aerodynamic  testing,  the  arcjet  flow  characteristics  need  to 
be  known  in  order  to  determine  if  they  provide  optimum  simulation  of  the  actual 
flow  conditions. 

In  arcjet  operations,  extremely  high  temperatures  are  generated.  For  example, 
for  the  test  reported  here,  the  stagnation  temperatures  are  of  the  order  of  6000*K 
(10,800°R).  Probes  must  be  highly  cooled  in  order  to  survive  the  environment. 
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Conventional  techniques  using  thermocouples,  for  example,  to  determine  local 
stream  energy  simply  will  not  work  because  of  these  high  temperatures.  Stream 
tube  flow-swallowing  techniques  can  be  used  up  to  moderate  temperatures, 
however,  such  enthalpy  probes  require  sharp  leading  edges;  and  these  can  not  be 
cooled  in  the  arcjet  environment.  Moreover,  such  probes  with  blunt  cooled  leading 
edges  yield  an  uncertain  stream-tube  capture;  and  resolution  and  accuracy  both  have 
large  uncertainties.  To  circumvent  these  problems  for  arcjet  stream  surveys,  a  simple 
method  has  been  developed  that  uses  small  diameter,  cooled,  hemispherical-shaped 
probes  to  measure  local  heat  transfer  and  pitot  pressure.  It  is  shown  in  this  report 
that  these  two  measurements  yield  stream  profiles  to  a  high  degree  of  resolution. 
The  theoretical  method  derived  here  has  been  formulated  into  a  computer  code  for 
data  reduction  of  stream  surveys.  This  code  can  be  incorporated  into  a  real-time  data 
acquisition  system  and  will  provide  on-line  reduced  stream  profiles  during  the 
nozzle  survey  for  immediate  assessment. 

PROFILE  is  a  computer  code  developed  to  use  the  60-MW IHF  stream  survey 
data  for  a  13-inch  conical  nozzle.  The  PROFILE  code  directly  produces  normalized 
velocity,  density,  mass  flux  and  enthalpy  profiles.  The  derivation  of  PROFILE  uses 
the  basic  conservation  equations  and  is  described  in  Chapter  IV.  PROFILE  uses  the 
heat  transfer  and  stagnation  pressure  data  as  inputs  and  outputs  normalized  values  of 
heat  transfer  rates,  stagnation  pressure,  velocity,  density,  mass  flux  and  enthalpy 
with  respect  to  centerline  values  of  the  flow.  Characteristics  of  flows  for  other  arcjet 
designs  and  other  nozzles  can  be  readily  obtained  using  heat  transfer  and  stagnation 
pressure  survey  data  in  conjunction  with  the  PROFILE  code  developed  here. 
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II.  INTERACTION  HEATING  FACILITY  DESCRIPTION 


A  schematic  diagram  and  photograph  of  the  IHF  arcjet  facility  are  shown  in 
Figure  1  and  2,  respectively.  The  major  components  of  the  facility  include  a 
constricted-arc  air  heater,  a  150  MW  power  supply,  two  interchangeable  hypersonic 
Mach  number  nozzles,  a  walk-in  test  chamber,  and  ancillary  subsystems  consisting 
of  a  steam-ejector  vacuum  system,  cooling  water  system,  and  automated  data 
acquisition  system.  The  high-power  capability  of  the  arc  heater  (rated  to  60  MW) 
provides  for  full  scale  subsystem  tests  in  high  enthalpy  streams  in  both  large  scale 
stagnation  flows  and  boundary  layer  flows  at  high  Reynolds  numbers. 

The  primary  component  of  the  facility  is  a  constricted-arc  heater  390  cm 
(153.3  in.)  long  and  8  cm  (3.15  in.)  in  diameter  shown  in  Figure  3.  The  rated 
airflow  rate  is  1 .4  kg/s  (3.09  lbs/s)  and  the  rated  pressure  is  12  atm.  Cooling  water 
is  supplied  at  7000  kPa  (1000  psi),  which  provides  cooling  for  a  maximum  heating 
of  12  kW/cm2  to  the  constrictor  wall.  Eight  anode  and  cathode  electrodes,  each 
capable  of  conducting  1000  A,  provide  a  maximum  arc  current  capacity  of  8000  A. 
Thus,  the  capability  for  continuous  operation  at  high  current  and  high  enthalpy 
levels  is  provided. 

Direct  current  for  operating  the  arc  heater  is  provided  by  a  150  MW  D.C. 
power  supply.  Six  modules  make  up  the  power  supply,  each  consisting  of  a  three- 
phase,  full-wave,  phase  controlled  silicon  rectifier.  One  module  is  capable  of 
operating  up  to  a  level  of  5500  V  into  an  arc  load  at  current  levels  up  to  2700  A. 
The  six  modules  can  be  connected  in  any  compatible  series/parallel  arrangement  by 
means  of  remotely  operated  setup  switches.  The  duty  cycle  for  the  power  supply  is 
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30  minutes  on,  30  minutes  off.  A  150  MW  rating  is  achieved  for  run  durations  of 
approximately  15  seconds. 

The  facility  has  two  interchangeable  nozzles:  1)  a  semielliptic,  Mach  number  5 
nozzle  with  an  open-jet,  flat-plate  test  section  measuring  80  X  80  cm  (31.5  X  31.5 
in.)  for  boundary-layer  surface  flow  tests  and  2)  a  conical,  nozzle  for  testing  large- 
scale  stagnation  flows  in  a  free  jet.  The  conical  nozzle  is  constructed  in  5  sections 
with  exit  diameters  varying  from  15.24  cm  (6  in.)  to  105  cm  (41.4  in.).  The  Mach 
number  varies  with  the  nozzle  area  ratio  from  3  to  7.5  depending  on  the  exit 
diameter  selected.  Both  nozzles  are  water  cooled  by  internal  water  passages,  and  the 
throat  diameters  are  6.03  cm  (2.37  in.).  Both  nozzles  are  designed  to  flow  into  an 
evacuated,  walk-in  test  chamber  that  can  be  maintained  at  a  pressure  of  less  than 
0.004  atm  (3  mmHg),  a  pressure  that  varies  with  flow  rate  according  to  the 
pumping  curve  of  the  steam  ejector  system.  Tests  have  shown  that  supersonic  flow 
with  33  %  blockage  ratio  has  been  attained  in  the  conical  nozzle,  corresponding  to  a 
model  diameter  of  about  60  cm  (24  in.). 

An  on-line  automated  data  acquisition  and  computing  system  is  used  to 
record  arcjet  performance  and  test  data.  The  conversion  from  analog  to  digital 
signals  is  made  by  means  of  an  integrating  digital  voltmeter,  so  that  the  signals  can  be 
processed  by  an  on-line  computer.  Output  consists  of  a  real-time  display  of  five 
channels  that  are  monitored  during  a  run,  printed  and  reduced  data,  and  a 
continuous  time  history  data  plot.  [Ref.  1] 
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m.  SURVEY  DESCRIPTION 


The  design  of  the  IHF  allows  for  the  control  of  two  variables:  1)  constrictor 
pressure  and  2)  arc  current.  The  test  conditions  are  uniquely  determined  by  these 
two  control  settings  only.  For  this  survey,  the  test  chamber  pressure  was 
maintained  at  a  low  pressure  to  allow  the  measurements  to  be  made  in  a  free  jet.  No 
attempt  was  made  to  control  the  expansion  ratio  into  the  free-jet  chamber.  In  general, 
the  flow  was  as  an  under-expanded  free  jet.  This  type  of  operation  causes  the  free  jet 
to  spread  and  to  operate  at  a  slightly  higher  Mach  number  relative  to  the  nozzle  exit 
plane.The  survey  measurements  close  to  the  nozzle  exit  plane,  then  are  relatively 
uninfluenced  by  the  free-jet  expansion  process.  Most  of  the  discussion  pertains  to 
the  7.62-cm  (3-inch)  station  to  avoid  the  complexity  of  the  additional  expansion. 
Constrictor  pressure  was  varied  for  several  values  of  current  in  order  to  map  the 
operating  envelope  for  stagnation  flows  using  the  conical  nozzle. 

The  probes  used  to  obtain  data  consisted  of  a  hemispherical  Gardon  type 
calorimeter  and  a  pitot  probe  shown  in  Figures  5  and  6.  The  probes  were  1 .58  cm 
(0.625  in)  in  diameter  and  were  mounted  20.32  cm  (8.0  in)  apart  horizontally.  The 
probes  were  mounted  to  a  transverse  mechanism  that  was  positioned  to  pass 
horizontally  through  the  flow.  Sweep  measurements  were  taken  at  7.62  cm  (3.0 
in),  19.5  cm  (7.7  in),  34.8  cm  (13.7  in)  and  52.6  cm  (20.7  in)  downstream  of  the 
nozzle  exit. 

Figure  7  shows  probes,  probe  supports  and  their  relationship  to  the  nozzle 
during  a  test  using  a  15.24  cm  (6.0  in.)  exit  diameter  nozzle.  Figure  8  shows 
probes  entering  the  flow.  Figure  9  shows  probes  approximately  on  centerline  of 
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nozzle.  All  photographs  were  taken  using  only  radiated  light  from  gas  behindthe 
normal  shocks  on  probes  and  support  mechanisms. 

Data  consisted  of  heat  transfer  rate  and  stagnation  pressure  measurements  taken 
at  0.212  cm  (0.083  in.)  intervals  across  the  flow.  Typical  profiles  for  heat  transfer 
rate  and  stagnation  pressure  for  each  current  setting  are  shown  in  Figure  10. 


IV.  DERIVATION  OF  DIRECT  METHOD  FOR  PROFILES 


Determining  the  stream  profiles  for  an  arcjet  wind  tunnel  nozzle  requires  that  the 
heat  transfer  rate  and  the  stagnation  pressure  measurements  behind  a  normal  shock 
make  up  the  only  parameters  on  which  the  velocity  characteristics  depend.  Briefly, 
this  method  is  based  on  the  fact  that  stagnation-point  heat  transfer  is  proportional  to 
the  product  pV3  (approximately)  while  the  pitot  pressure  is  proportional  to  pV2. 
When  these  two  quantities  are  known,  a  simultaneous  solution  gives  density  and 
velocity  directly  from  the  measurements.  Once  the  velocity  profile  has  been 
determined  the  profiles  for  enthalpy,  mass  flux  and  density  can  be  readily  found. 

In  the  following  derivation,  free  stream  flow  parameters  are  those  found  in  front 
of  the  normal  shock  and  will  be  designated  by  subscript  one  while  subscript  two 
will  denote  parameters  behind  the  normal  shock.  The  starting  point  is  the 
momentum  and  continuity  equations  in  the  following  form: 

+  a) 

PlVl=PjV2  (2) 

Rearranging  the  momentum  equation  and  substituting  into  it  the  continuity  equation 


gives: 


=p^‘-p?r 


Solving  for  p2 


P2~P"  l«-(i Wp,)) 
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Noting  that  M2  is  much  less  than  1 ,  that  M2  =  <p(  Y»  M)  and  recalling  that  from  state  2 
to  stagnation  state  2  is  along  an  isentrope  then: 


Normalizing  equation  (5)  with  centerline  values  results  in 


Pi 
( Pi ) 


Substituting  equation  (6)  into  a  normalized  equation  (4)  gives: 


(pv>. 


p,v, 

(p.v?) 


At  this  point  the  normalized  stagnation  pressure  behind  the  shock  is  related  to 
the  velocity  in  front  of  the  shock.  The  next  step  introduces  the  closed  form  of  the 
Fay-Riddell  [Ref  2]  equation  derived  in  Appendix  A.,  which  is  a  correlation  for  the 
heat  transfer  rate  to  a  hemisphere;  namely. 


9  =  0  J0417  /  — -  (ht-  A* )  - 

V  R  K  ‘  (Btu/sec-ft2)  (8) 


where  ht  is  the  total  enthalpy  of  the  flow  and  hw  is  the  enthalpy  of  the  wall.  The 
units  are:  h,  Btu/lbm;  pe,  atm;  and  R,  ft.  Normalizing  with  centerline  values  gives: 
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(9) 


(«) 


CL 


(P‘>) 


CL 


(  (*,-  M  \ 
(*,-  k)cl 


Using  the  following  form  of  the  energy  equation 


*i  + 


(10) 


and  equation  (7)  in  equation  (9)  while  neglecting  the  enthalpy  of  the  wall  compared 
to  the  total  enthalpy  of  the  flow  yields: 


4  I  p  ,v,*  (y,‘+z*,) 

<«)„'  J(rrf)'L  O'.1 + 

Rearranging 


(ID 
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(4) 


CL 


( 

1  + 
v _ 

( 

1  + 

V 


Uh) 

v 2 
vl  l 


2\\ 


i  J 


CL 


(12) 


Since  the  velocity  in  front  of  the  shock  is  supersonic  the  last  term  drops  out. 
Changing  notation  to  reflect  only  centerline  and  free  stream  gives: 
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(13) 


Squaring  equation  (13)  and  dividing  by  equation  (7)  removes  the  density  term  and 
results  in  an  equation  which  relates  free  stream  velocity,  stagnation  pressure  and  heat 
flux. 


(14) 


This  equation  equates  the  normalized  velocity  of  the  free  stream  to  the  two  measured 
quantities  behind  the  normal  shock. 

The  assumption  made  to  derive  this  result  are  simply  that  the  flow  is 
axisymmetric  and  the  enthalpy  of  the  wall  is  small  compared  with  the  total  enthalpy 
of  the  free  flow.  The  wall  enthalpy  is  measured  directly  from  cooling  water  and  is 
found  to  be  approximately  0.3  MJ/kg  (134  Btu/lbm)  or  less  than  5  %  of  total 
enthalpy. 

Once  the  normalized  velocity  has  been  found  it  is  used  in  equation  (7)  to  find  the 
value  for  the  normalized  density. 


(15) 


Normalized  mass  flux  is  obtained  from  the  product  of  normalized  density  and 
velocity. 


juJL 

ntf. 


CL 


(16) 


The  normalized  enthalpy  is  taken  directly  from  the  energy  equation  where 
gravitational  effects  are  neglected  and  the  static  enthalpy  is  assumed  to  be  much  less 
than  total  enthalpy,  which  results  in  the  following  equation. 


10 


CL 


CL 


(17) 

Equations  (14),  (IS),  (16),  and  (17)  were  used  in  the  code  PROFILE  to  calculate  the 
characteristic  profiles  of  the  arcjet  high  enthalpy  flows  using  survey  data. 


V.  RESULTS 


For  high-powered,  arc-heated,  hypersonic  facilities  the  testing  for  aerodynamic 
or  thermodynamic  (heat  transfer)  studies  requires  a  central  core  flow  with  properties 
well  characterized.  Ideally,  of  course,  constant  core  flow  properties  are  desirable. 
The  present  study  was  undertaken  to  characterize  the  flow  properties  of  the  60- MW 
arc  heated  wind  tunnel  at  NAS  A- Ames  Research  Center.  Such  flow  properties  are 
essential  to  the  engineer  for  interpretation  of  the  results  for  aero-thermodynamic 
tests. 

Arc  heaters  operate  with  a  large  enthalpy  gradients  near  the  walls.  These 
gradients  cause  the  flow  properties  in  the  central  core  to  vary;  and  are  an  important 
concern  in  simulation  in  these  ground-based  test  facilities.  The  degree  of  this 
variation  has  not  been  reported  up  to  the  present  time  for  the  arcjet  facility  at  NASA- 
Ames  Research  Center.  Enthalpy  gradients  are  inherent  in  all  types  of  arc  heaters 
that  operate  with  cooled  walls.  For  example,  the  temperature  ratio  of  the  wall  for  a 
typical  arc  heater  varies  from  2  -  5  %  when  the  enthalpy  level  is  in  the  typical  vehicle 
reentry  range  of  18.6  -  41.87  MJ/kg  (8000-18000  Btu/lbm)  which  in  terms  of 
velocity  equates  to  20,000-30,000  ft/sec.  The  constricted  arc  heater  was  developed 
to  minimize  this  gradient  by  the  current  density  and  length  to  cause  the  enthalpy 
profile  to  assume  a  flatter  shape  than  a  free-burning  or  vortex-stabilized  arc  heater. 
The  theory  of  Stine-Watson  [Ref.  3]  shows  that  for  a  constricted  arc  heater  the 
enthalpy  profile  assumes  a  Bessel  function  form  as  opposed  to  a  linear  profile  for  a 
free-burning  or  vortex-stabilized  design  (Huels  arc)  illustrated  in  Figure  1 1 . 

Results  from  the  profile  measurements  obtained  here  can  be  used  to  answer  the 
question  as  to  whether  the  property  variables  (enthalpy,  velocity,  density  and 
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energy  flux)  have  the  same  profile  or  shape  at  the  nozzle  exit  as  they  do  for  the 
constrictor.  Scant  measurements  exist  for  the  constrictor.  Most  of  the  values  for  the 
constrictor  arc  heaters  are  predicted  by  the  Sdne- Watson  theory.  Elaboration  of  this 
theory  is  represented  by  a  computer  code  called  ARCFLO  II  [Ref.  4].  This  code 
was  developed  to  predict  the  performance  of  high-pressure,  high-enthalpy 
constricted-arc  heaters.  The  conservation  equations  (mass,  axial  momentum,  and 
energy)  are  written  for  a  steady-state  and  radiatively  participating  gas  mixture  with 
electric  current  conduction.  The  governing  equations  employ  the  boundary-layer 
approximation  and,  hence,  the  pressure  in  the  wall  normal  direction  is  assumed 
constant.  Auxiliary  relations,  such  as  the  equation  of  state,  global  current  continuity, 
and  global  mass  continuity,  are  used  to  relate  the  gas  density  with  the  pressure,  mass 
flow  rate,  and  mixture  properties. 

The  profiles  of  density,  velocity,  enthalpy  and  energy  flux  computed  by 
ARCFLO  II  in  the  constrictor  are  compared  with  experimental  survey  data  obtained 
at  the  nozzle  exit  in  Figure  12.  The  difference  seen  in  the  density  profiles,  Figure 
12a,  can  be  accounted  for  by  the  increase  in  density  near  the  wall  of  the  constrictor  as 
opposed  to  the  decrease  in  density  in  an  expanding  free  jet.  The  velocity  profiles, 
Figure  12b,  are  very  similar  with  the  experimental  data  showing  a  slightly  flatter 
profile.  The  enthalpy  profiles,  Figure  12c,  shows  for  ARCFLO  II  a  central  core 
region  where  the  enthalpy  drops  moderately  from  the  centerline  to  r/R  =  0.4,  a 
linearly  decreasing  region  to  r/R  =  0.95,  and  then  a  steep  drop  near  the  wall  where 
the  conduction  losses  dominate.  The  ratio  of  centerline  enthalpy  to  mass  average 
enthalpy  is  1 .6.  This  low  ratio  value  is  typical  for  wall-stabilized  arcs.  The  survey 
data  exhibits  a  much  flatter  profile  than  ARCFLO  II  with  a  ratio  of  centerline 
enthalpy  to  mass  average  enthalpy  of  1.15. 
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The  product  of  the  above  three  profiles  is  shown  in  Figure  12d  and  is  the 
energy  flux  of  the  airstream.  Note  that  this  parameter  is  proportional  to  the 
convective  heat-transfer  rate  to  a  body  within  the  stream  and  that  it  is  a  measure  of  die 
usable  stream  area.  The  flat  profile  for  ARCFLO  n  implies  that  stagnation  point 
heating  to  a  test  body  would  be  constant  to  within  10  %  of  the  value  obtained  on  the 
centerline  flow  for  r/R  up  to  0.7.  The  survey  data  profile  shows  a  more 
pronounced  curve  which  can  be  attributed  to  die  decrease  in  density  in  the  free  jet.  If 
test  chamber  static  pressure  is  maintained  at  the  nozzle  exit  static  pressure  then  the 
density  profile  would  become  flatter  and  thereby  improve  the  energy  flux  profile. 

Figures  13  -  28  are  plots  of  PROFILE  code  output  data  at  specified  survey  test 
points.  All  test  point  plots  are  normalized  with  respect  to  the  corresponding  values 
obtained  at  the  3-inch  downstream  positions. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 

The  present  study  of  survey  data  for  the  13-inch  IHF  nozzle  showed: 

That  there  is  good  agreement  between  the  ARCFLO II  code 
prediction  for  velocity  profile  and  experimental  data; 

Enthalpy  profile  of  experimental  data  is  much  flatter  than  ARCFLO  II 
theory  predicts,  which  can  be  explained  by  the  codes  inability  to 
account  for  energy  transfer  toward  constrictor  edges. 

It  is  recommended  that: 

Future  surveys  control  test  chamber  conditions  in  order  to  determine 
if  improvements  in  the  profile  characteristics  result; 

Sting  mounted  heat  transfer  and  stagnation  pressure  probes  be 
developed  in  order  to  reduce  probe  mount  interference  with 
boundary  layer  of  the  flow; 

Survey  data  be  organized  into  accessible  files  so  that  characteristic 
maps  for  individual  nozzles  and  complete  nozzles  can  be  utilized  by 
researchers  and  contractors. 
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APPENDIX  A 

DERIVATION  OF  THE  CLOSED  FORM  OF  THE  FAY-RIDDELL  EQUATION 

by  WARREN  WINOVICH 


1.  The  Nusselt  number  is  defined  by: 


|| 

Nu=-jj-L  (w  =  wall  value) 


(Al) 


2.  The  heating  rate  is  given  by: 


q  =  H  (T ,  -Tw) 


H(h.-hw) 


c'.c 


where 


(A2) 


c  = 


= ml 
V  _  m  m  *■ 


p  T  —  T 

r9  1  I  1  a 

(by  definition  of  Nusselt  number;  wall  temperature  used  as  reference) 

3.  Solving  for  heating  rate  in  terms  of  Nusselt  number. 


Nu  kw  ( h ,  -  hm ) 
9  ~  x  c,  c 


q  = 


m  c 


where  subscript  "e"  indicates  external  flow  conditions 


(A3a) 


(A3b) 


6 


Nu  rr-r. —  (*»  —  *») 


(A3e) 


4.  For  the  laminar  boundary  layer  near  the  stagnation  point: 


therefore 


(A4a) 


9  = 


Nu 

\fRe 


(»,-  «o 

*A"’  Pr.  c 


(A4b) 


5.  Evaluation  of  velocity  gradient  at  the  stagnation  point  for  a  hemisphere: 


pt  +  y -uj  =  pt  (momentum  equation  ) 


m 


(A5a) 


P  t  /  J f  \ 

-jp-  =  l-sin2(^J  (Newtonian  flow  over  hemisphere) 


(A5b) 


(A5c) 
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(4H-±,RS)‘ 


P.  =  P,  for  x  =  0 


(A5d) 


6.  The  heating  rate  equation  becomes: 


Vp.p. 


1 

VF 


( h ,  -  *,) 

Prm  c 


1.  The  Fay  and  Riddell  result  for  the  heating  transfer  coefficient  is: 


Nu 


=  0.763 (Pr) 


Ip.p-J  1 


Vr7 

(equation  63  and  45  Fay  and  Riddell) 

8.  Substituting  equation  (A7)  into  equation(A6);  using 


rv 
\p,<  > 


=  T,Z,  and  R  =  535 


gives: 


q  =  0.763  (/*r)  '  ^2gc(53. 3) 

Ip-pJ  V7*Z‘  V* 


rSjL  =  0.763 

-  h*>)  ~  (pr)~*  * 


x».l 

^2gc(53.3)VpX(^f)  </TX 


(A6) 


(A7) 


(A8a) 


(A8b) 
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9.  Using: 


P* 


39.65  tJ 


(A9a) 


P. 


=184x 10 


uoooJ 


•  .<21 


K. 


(A9b) 


where 


K  = 


Ptuct 

P  ferai  U 


Equation  (A9b)  is  appropriate  for  temperatures  of  6000*K  (10000’R)  where  arcjet 
operate  and  for  reentry  and  reentry  stagnation  conditions  (Vcntry  >  5000  m/sec  ~ 
15000  ft/sec) 


P.P.  = 


3 9.6 5 x  184 x  10" 7  PJ, 


•.<21 


1000 


«.<2t 


T.Z, 


K, 


P»P « 
P. 


«.<2( 

=  9.60  x  10*  -=f-=~  K, 


(A9c) 


„  „  „  *.<2* 

PePe  =  T,Z,  7\ _ 

p.p,  r.z,  r*«*  k. 

Note  that  Kw  =1.00  for  the  expression  chosen  for  viscosity. 

10.  The  final  form  for  the  expression  for  heating  rate  becomes: 


(A9d) 


gVR 


0.763 


•”  *.314 


Vp7(h.~  *•)  (/>r) 


rr«y2gc(53.3)  3.10x  10 


AT.z.)  t; 


( T,Z ,) 


•.5 


(AlOa) 
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Simplifying  for  Pr  =  0.71  the  equation  becomes: 


«v*  .  .  - 

—7=  /. - TT  =  (constant  } - 

*m)  (TmZm) 


(t.z.)"st:isik:4 

•.•372 


qy/~R 

Vp7(h,~  hw) 


=  (constant  } 


(T.Z.)  fZ.\ 

,  •.•372  I  7  I 

(rwzw) 


Noting  that  Zw  =  1.00 

For  Pr  =  0.71  the  constant  in  brackets  above  becomes: 


(constant  }  =  JL16jL  ^2gc(53.3)  3.10x  10'3  =  0.02225 

(Pr) 

and  the  expression  for  heating  rate  is: 


gV* 

^/~p7(h,~  hK) 


{0.02225}  \K 


*  \TmZj 


•  ■•372 


(A10b) 


The  term  within  the  braces  is  only  a  function  of  the  stagnation  state.  For  usual  wind 
tunnel  tests  the  stagnation  state  is  within  the  limits: 

0.01  <  ps  <  1  atmosphere  130  <  hs  <  20,000  Btu/lbm 
For  this  range  of  pressure  and  enthalpy  the  term  in  the  braces  is  essentially  a 
constant: 


=  1.873±  0.042 
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Therefore  the  hemisphere  formula  for  heating  rate  is: 


or: 


hw) 


{0.02225K1 .873} 


0.0417 


q  =°.«417^/A(*,-  A.) 

where 

q  =  Btu/s  ft2 

p,  =  atm  pt  (measured) 

R  =  feet 
h  =  Btu  /lb  m 


(AlOc) 


Symbols 

Cp  Specific  heat 

h  Entahlpy 

k  Thermal  Conductivity 

Pr  Prandtl  Number 

R  Radius 

T  Temperature 

Z  Compressibilty  Factor 

p  Density 

p  Viscosity 


Subscripts 

e  External  Conditions 

s  Stream  Conditions 

t  Stagnation  Conditions 

w  Wall  Conditions 
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APPENDIX  B 


UNCERTAINTY  ANALYSIS  FOR  DERIVATIONS 

Uncertainty  analysis  of  the  heat  transfer  and  stagnation  pressure  probes  is  based 
on  manufacture's  claims  for  uncertainties.  In  the  case  of  the  heat  transfer  probe  the 
manufacturer  list's  accuracy  of  readings  between  3  &  4  %  of  the  value  recorded. 
For  the  following  analysis  the  value  of  5  %  is  used.  The  accuracy  of  the  stagnation 
pressure  probe  has  been  found  by  NASA  personnel  to  be  on  the  order  of  1/4  %  of 
the  recorded  value.  Using  these  values  the  uncertainty  of  the  normalized  velocity  is 
shown  below. 

Starting  with  equation  (14): 


and  taking  the  logarithmic  differential  gives: 


Integrating 

d  V  1  f  d  4  d  icL  _  1  f±P  _  dPci) 

VCL  "  2  l  *  "  iCL  )  '  4  ‘  PCL  J 

v  =  ~2  [(e  r)^  “  (e  r)9ct]  "J  [(*0^  —  (®  r)  Pci] 

CL 

where  (er)  is  uncertainty. 

In  the  worst  case  where  centerline  values  are  opposite  in  sign  of  the  original 
values  results  in 
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=  y  [(0.05) -(-0.05)]  -  y  [(0.0025)  -(-0.002$)] 

CL* 


P~  =  4-  to  JO]  -4-  10.0051 


so  that 


±  5  % 


Applying  the  velocity  uncertainty  to  the  other  property  equations  (15),  (16)  and 
(17)  gives: 

« -«)■*••• 


■  mKfti 


±10% 
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APPENDIX  C 


COMPUTER  PROGRAM  FOR  DETERMINING  PROFILES 

FILE*  PR0F2  FORTRAN  A1 


PROGRAM  PROFILE 

THIS  PROGRAM  IS  DESIGNED  TO  CALCULATE  FLOH  PROFILES  FOR 
NASA  AMES  IHF  ARCJET  MIND  TUNNELS.  VELOCITY,  DENSITY,  MASS  FLUX 
AND  ENTHALPY  PROFILES  NORMALIZED  TO  CENTERLINE  CONDITIONS 
ARE  DEVELOPED. 


DIMENSION  QDOT(IOOO),  QDOTN(IOOO),  PTOTUOOO),  VNOR(IOOO), 

8  RHOUOOO),  MFLUX(IOOO),  HTOT(IOOO),  DIST(IOOO), 

8  QSQR(IOOO),  PTOTNC 1000) 

REAL  VNOR, RHO, MFLUX, HTOT,QSQR,PTOTN,QDOTN 
INTEGER  S,I, J,K,L, IMAX, IQDOT, I PTOT, VAL, 11,12 
REAL  PMAX,  QMAX, VMIN 
CHARACTER  *36  G 
CHARACTER  *1  E,Y,N 
CHARACTER  *6  A,B,C 
CHARACTER  *A5  F 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 
THIS  SECTION  COPIES  TEST  CONDITIONS  INTO  OUTPUT  FILE 
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXKXXXXXXXXXXKXXXXXXXKXXXXX 
READ(5, 1 )F 
FORMAT(AA5) 

WRITEC7,1)F 
READ(5,2)G,C 
FORMAT(A36,A6) 

MRITE(7,*)G,C 
DO  10  1=1,3 
READ(5,1)F 
WRITE( x, 1  )F 
WRITE(7,1)F 
0  CONTINUE 

XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXXKXXXXXKXXXXXXXX 
THIS  NEXT  SECTION  READS  TEST  DATA  FOR  QDOT  AND  PTOT. 

FIRST.  NORMALIZES  VALUES  OF  QDOT  AND  PTOT 
SECOND.  LOCATES  MAXIMUM  VALUE  TO  USE  AS  CENTERLINE 
THIRD.  CALCULATES  VALUES  FOR  VELOCITY,  DENSITY,  MASS  FLUX  AND 
ENTHALPY  PROFILES  NORMALIZED  TO  CENTERLINE 

NOTE.  QDOT  AND  PTOT  DATA  SHOULD  CONTAIN  ONLY  ONE  ZERO  AS  LAST 

VALUE.  THIS  IS  CRITERIA  USED  TO  TERMINATE  READING  OF  DATA. 

xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

IMAX=0 

DO  20  1=1,1000 

READ(5,15)QD0T(I),  PTOT(I),  DIST(I) 

IF( QDOT ( I ) . EQ . 0 )GOTO  21 
IMAX=IMAX+1 

20  CONTINUE 

21  WRITE!*,*)’ NUMBER  OF  DATA  POINTS  READ  =’,IMAX 
ENDFILE! 5) 
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C  XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXK 

C  1  NORMALIZATION  OF  QDOT  AND  PTOT 

C  XXXXXXXXXXXXXXXXXXXXXXXXXKXXXXXXKXXKXXXXXXXXKXKXKXXKKXXXXKXKKKX 

QMAX  =  0 
PMAX  =  0 

DO  200  I  =  1 , I MAX 

QMAX  =  MAX(QMAX,QDOT(I>) 

PMAX  =  MAX( PMAX, PTOT ( I ) ) 

200  CONTINUE 

C  XXXXXXXXXXXXXXXKXXXXXXXXKXXXXXXXXXXXXXXXXXXXXXXKKXXXXXKXXKXXXXXX 

C  2  DETERMINING  LOCATION  OF  MAX  VALUES.  THESE  HILL  BE  USED  TO 

C  ALIGN  PROBE  VALUES  TO  CENTERLINE. 

C  XXXXXXXXXXXXXXXXXXXXXXXXXXXXKXXXXXXXXXXXXKXXXXXXXKXXXKXXKXXKXXKXX 

DO  201  1  =  1 , IMAX 
IQDOT=I 

VAL=DIMCQMAX,QDOT(I)) 

IFCVAL . EQ . OIGOTO  202 

201  CONTINUE 

202  DO  203  1=1, IMAX 

IPTOT= I 

VAL 1 =DIM(PMAX, PTOT ( I ) ) 

IFCVALl .EQ.OIGOTO  205 

203  CONTINUE 

205  HRITEC*,#)  'I  AT  QMAX  =',IQD0T,*  I  AT  PMAX  *MPTOT 

C  XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

C  3  CALCULATION  OF  PROFILES 

C  xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

Il=IQDOT+179 
I2=IPTOT+179 

WRITECX, X) 'll  II, '  12  =',I2 

K  =  0 

DO  300  I=IQDOT , II 
K  =  K+1 

QDOTN(K)=ABS(QDOT(l)/QMAX) 

QSQR(K)=QD0TN(K)XX2 

300  DISTCK)=DISTCI 1-DIST  L I QDOT ) 

K  =  0 

DO  301  I=IPT0T , 12 
K  =  K+1 

301  PTOTN(K)=ABS(PTOTCI )/PMAX) 

DO  302  1=1, K 

VNOR(I)  =  SQRT(SQRT(QSQRCI)/PTOTNCI))) 

RHOCI)  =  PT0Tmi)/U/VN0R(l))XK2 
MFLUX(I)  =  RHOmxVNORII) 

302  HTOT(I)  =  VN0R(I)XX2 
GOTO  A00 


A00  HRITE(7,X)'QMAX  =  *,QMAX 
WRITEC6 ,  X ) ' QMAX  =  * , QMAX 
HRITE(7,X) 'PMAX  =',PMAX 
WRITEC  6 , X) ' PMAX  =’,PMAX 

HR I T  E ( 7 , X ) *  a  s  s  s  a  a  s  a  s  s  s  a  r  a  a  a  a  a  a  s  a  a  a  a  a  a  a  a  a  s  a  a  s  s  a  a  a  a  a  a  s  a  a  s  a  a  s  s  a  s  a  a  a  s  r 

B================* 

WRITEC  7 , X) '  HEATING  STAG  VELOCITY  DENSITY  MFLUX 
6ENTHPY  DIST ' 

WRITE( 7 , * ) '  RATE  PRESSURE* 

WRITE( 7,X)'=r5sarrars8aBsarrsr=sa=sa=saas=aaraas8asaaase88srsaasss 
ftaasaaaaaaaaaaaaa* 

MRITE(7,600)( QDOTN( I ) » PTOTNC I ) , VNOR( I ) , RH0( I ) ,MFLUX( I ) , HTOT( I ) , DIS 
BTC  I ) , 1  =  1 , K, A) 

600  FORMAT(6(5X(F5.3)),5X,F6.3) 

STOP 

END 
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OPERATOR 
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Schematic  diagram  of  the  60  MW  IHF, 


Photograph  of  the  60  MW  IHF  with  Mach  7.5  conical  nozzle. 


ARGON  START/SHIELD  GAS 
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CONSTRICTED  ARC  HEATER 
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Figure  5  Gardon  type  calorimeter  probe. 


Figure  6  Stagnation  pressure  probe  -  with  support  mount. 
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Figure  7  Probes  and  support  mounts  outside  6-inch  nozzle  free  jet  flow. 


Figure  8  Probes  entering  6-incli  nozzle  free  jet  Flow  during  survey. 
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Figure  9  Probes  near  centerline  of  6-inch  nozzle  free  jet  flow. 


1600  Amps 
2000  Amps 
2500  Amps 
3000  Amps 
4000  Amps 
5000  Amps 
6000  Amps 


(a)  Heat  transfer  rate  verses  radius  at  3-inch 
position  from  nozzle  exit 


1600  Amps 
2000  Amps 
2500  Amps 
3000  Amps 
4000  Amps 
5000  Amps 
6000  Amps 


(b)  Stagnation  Pressure  verses  radius  at  3-inch 
position  from  nozzle  exit 


Figure  10  Typical  profiles  of  the  heat  transfer  rate 
and  stagnation  pressure,  13-inch  nozzle  (A/A*  =  30) 
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Figure  11  Comparison  of  enthalpy  profiles  for  8-cm 
constricted  arc. 
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Figure  12  Comparison  of  ARCFLO  II  and  survey  data  profiles  of 
8-cm  constricted  arc. 
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PROFILE 


Figure  15  13-inch  nozzle  profiles  normalized  by  3-inch  position  from  exit. 

Constrictor  conditions:  2000  Amps,  46  psia  constrictor  pressure,  7.4  MW  arc  power  input. 


PROFILE 
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Figure  16  13-inch  nozzle  profiles  normalized  by  3-inch  position  from  exit. 

Constrictor  conditions:  2500  Amps,  28  psia  constrictor  pressure,  9 3  MW  arc  power  input. 


VELOCITY  PROFILE  '  L _  ENTHALPY  PROFILE 
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Figure  18  13-inch  nozzle  profiles  normalized  by  3-inch  position  from  exit. 

Constrictor  conditions:  2500  Amps,  72  psia  constrictor  pressure,  93  MW  arc  power  input. 
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Figure  22  13-inch  nozzle  profiles  normalized  by  3-inch  position  from  exit. 

Constrictor  conditions:  4000  Amps,  30  psia  constrictor  pressure,  14.8  MW  arc  power  input. 
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Figure  23  13-inch  nozzle  profiles  normalized  by  3-inch  position  from  exit. 

Constrictor  conditions:  4000  Amps,  51  psia  constrictor  pressure,  14.8  MW  arc  power  input. 
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Figure  26  13-inch  nozzle  profiles  normalized  by  3-inch  position  from  exit. 

Constrictor  conditions:  5000  Amps,  53  psia  constrictor  pressure,  18.5  MW  arc  power  input 
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Figure  27  13-inch  nozzle  profiles  normalized  by  3-inch  position  from  exit. 

Constrictor  conditions:  5000  Amps,  78  psia  constrictor  pressure,  18.5  MW  arc  power  input, 
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Figure  28  13-inch  nozzle  profiles  normalized  by  3-inch  position  from  exit. 

Constrictor  conditions:  6000  Amps,  109  psia  constrictor  pressure,  22  MW  arc  power  input. 


LIST  OF  REFERENCES 


1.  Winovich,  W.,  Balboni,  J.  Balakrishnan,  A..,  "Application  of  Numerical 

Simulation  to  Enhance  Arcjet  Performance,"  AIAA  Paper  85-1006,  1985, 
pp.  398-400. 

2.  Fay,  J.  A.  and  Riddell,  F.  R.,  "Theory  of  Stagnation  Point  Heat  Transfer  in 

Dissociated  Air,"  Journal  of  the  Aeronautical  Sciences  ,  Vol.  25, 1958,  pp. 
73-85. 

3.  Stine,  H.  A.„Watson,V.  R.,  "The  Theoretical  Enthalpy  Distribution  of  Air  in 

Steady  Flow  Along  the  Axis  of  a  Direct-Current  Electric  Arc,"  NASA  TN 
D-1331,1962. 

4.  Nicolet,  W.  E.  and  others,  "Analytical  and  Design  Study  for  a  High-Pressure, 

High-Enthalpy  Constricted  Arc  Heater,"  AEDC  TR-75-47,  1975. 


INITIAL  DISTRIBUTION  LIST 


No.  Copies 


1 .  Defense  Technical  Information  Cento'  2 

Cameron  Station 

Alexandria,  Virginia  22304-6145 

2.  Library,  Code  0142  2 

Naval  Postgraduate  School 

Monterey,  California  93943-  5002 

3.  Chairman,  Code  67  1 

Department  of  Aeronautics  and  Astronautics 

Naval  Postgraduate  School 
Monterey,  California  93943-5360 

4.  Prof.  Richard  Wood,  Code  67Wr  2 

Naval  Postgraduate  School 

Monterey,  California  93943-5360 

5.  Prof.  Raymond  Shreeve,  Code  67Sf  1 

Naval  Postgraduate  School 

Monterey,  California  93943-5360 

6.  Warren  Winovich  Code  RTF  2 

NASA  Ames  Research  Center 

Moffett  Field,  California  94035 

7.  John  Balboni  Code  RTF  2 

NASA  Ames  Research  Center 

Moffett  Field,  California  94035 

8.  Robert  W.  Kopp  3 

National  Aero  Space  Plane 

AFSC/NAG 

Wright-Patterson  AFB,  Ohio  45433-6503 


55 


9.  Howard  Goldstein  CodeR 
NASA  Ames  Research  Center 
Moffett  Field,  California  94035 


10.  Jim  Arnold  Code  R 

NASA  Ames  Research  Center 
Moffett  Field,  California  94035 


11.  A1  Covington  CodeR 
NASA  Ames  Research  Center 
Moffett  Field,  California  94035 

12.  Tom  Polek  Code  RTF 
NASA  Ames  Research  Center 
Moffett  Field,  California  94035 


56 


